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SPECIFICATION 

TITLE OF THE INVENTION 

Methods of Designing Optimal Linear Controllers 

TECHNICAL FIELD OF THE INVENTION 

This invention relates to the design of optimal multivariable linear controllers and PID 
controllers. 

BACKGROUND OF THE INVNTION 

A PID (proportional-integral-derivative) controller is a device, usually implemented in a 
computer, that is used to control a process variable (e.g., temperature, pressure, etc) of 
an industrial process. For a multiple-input and multiple-output (MIMO) process, the 
process variable y(t) and controller output u(t) are n-dimensional vector variables of the 
time t. The controller receives the sample values of y(t) and calculates its output u(t) 
according to a PID equation. The controller sends its output u(t) to the process so that the 
error e(t)=r(t)-y{t) approaches zero as time t increases, where r(t) is the set point (also 
known as the reference signal or the command signal, etc). 

There are many types of PID controllers, depending on the use of different types of PID 
equations. In discrete time form, all these PID controllers can be viewed as special cases 
of the general linear controller with the linear control equation 

DUk=ErrCyk 

where yk=y(tk) is the process variable at time tk=to+kTs. to is the initial time, Ts>0 is the 
constant sampling period, k=0, 1, 2, ... is a non-negative integer called discrete time 
variable, Uk=u(tk) is the controller output at time tk, and Uk can be subject to upper limit 
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and/or lower limit constraints placed on one or more of its components, ri(=r(tk) is the set 
point at time tk, D, E and C are n by n matrix polynomials in the unit backward shifting 
operator z'^ such that for any discrete time signal Xk, z'^Xk=Xk.i, and one or more of the 
said D. E and C contain tuning parameters that are to be determined. 

For example, if D=(1- z"^)/Ts*l. where I is the identity matrix of order n. and 
E=C=Ki+KpD+KdD^, where Kj, Kp and Kd are user specified constant coefficient matrices 
called the integral gain, the proportional gain and the derivative gain, then the linear 
controller becomes the following "type A" PID controller 



where ek=rk-yk is the error. The three gains Kj, Kp and Kd are the tuning parameters. 

If the aforesaid E is changed to K|+KpD or Kj, then the linear controller becomes the 
following "type B" and "type C" PID controllers, respectively: 



Once its structure is selected, the performance of a linear controller or a PID controller 
depends mainly on the choice of its tuning parameters. 

How to properly choose the tuning parameters for a PID controller is a problem that has 
attracted a lot of studies ever since PID controllers became widely used in industry in the 
early 1940s. The Ziegler-Nichols tuning methods developed by Ziegler and Nichols in 
1942 and 1943 are still widely used in industry. Other model based methods choose the 
tuning parameters by minimizing some well-known control performance index such as the 



Duk= Kiek+ KpDek+ KdD^ek 



DUk=Kiek+KpDek-KdDVk 
DUk=Kiek-KpDyk-KdD^yk 
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integral absolute errors (lAE), the integral squared errors (ISE), etc. However, practice 
shows that all these methods often lead to oscillatory control results. 

DETAILED DESCRIPTION OF THE INVENTION 

Suppose the open-loop transfer function of the process from Uk to yi< is A"^B, where A and 
B are known n by n matrix polynomials in the unit backward shifting operator z'V 

A=l-AiZ"^-A2Z'^- ... -ApZ-P 
B= BiZ^ + B2Z^+ ... -HBqZ-^ 

where I is the identity matrix of order n, p and q are two known positive numbers, and Ai, 
A2, Ap, Bi, B2. .... and Bq are constant coefficient matrices of order n. A and B can be 
obtained by use of the commercially available software product "System Identification 
Toolbox" from The MathWorks Inc. Below is another way for finding A and B. 

Select a positive integer v that is not less than n(p+q) and let k be not less than v + max{p, 
q}, where max{p. q} denotes the larger one of the p and q. Let hk be a vector defined by 

hK'^=[yk./ yk.2^... yk/ Uk./ Uk.2'^... UkV I] 
where for any matrix X. denotes the transpose of X, and let 

Y=[yk yk-i ... yk.v]^ 
H=[hk hk-i ... hk-vr 
W=diag(1, w, w^, .... w") 

where w is a positive real number not larger than one, and W is a diagonal matrix with 1, 
w. w^ and w"" on its diagonal. The number k should be large enough and the controller 
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output signal Uk should be chosen so as to ensure that the matrix H is of full column rank. 
Then the matrices Ai, A2, ... , Ap, 61,62 and eq, and the vector d can be obtained from 

[Ai A2...Ap 61 82... ec, d] = ((WHr(WY))^ 

where for any matrix X, denotes the Moore-Penrose pseudo-inverse of X, which is a 
standard function under the name "pinv" in the software product MATLAe from The 
MathWorks Inc. (MATLAB is a registered trademark of The MathWorks Inc.). Thus the 
two polynomials A and 6 and a vector d are obtained. When A. 8, and d are so obtained, 
the process can be approximately represented by the following linear model: 

Ayk=Buk+d 

And the open-loop transfer function of the process is then A"^B. Please note that this 
invention does not deal with the problem of how to get the polynomials A and B. From 
now on it is always assumed that, one way or another, the two polynomials A and 6 are 
known. 

If the process is controlled by the following linear controller as described in the previous 
section: 

DuK=Erk-Cyk 

where D, E, and C are n by n polynomials in the backward shifting operator z \ and one 
or more of D, E and C contain tuning parameters that are to be specified, then the 
closed-loop transfer function from the set point rk to the process variable yk is 



Q=(A+BD-^C)'^BD-^E 
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From Q the characteristic polynomial in z, denoted by b(z), can be found. All roots of the 
characteristic equation b(z)=0. denoted by Zi, Z2, Z3. .... Zj. where j is a positive number, 
form all the poles of Q. 

Let max {|zi|. IZ2I, ... , |Z||} denote the maximum of all the absolute values |zi|, |Z2| and 

|Zj| of the q poles of Q. Then the best choice of the tuning parameters is such that the max 
{|zi|, IZ2I, .... |Zj|} is minimized, i.e., the best tuning parameters form a solution to the 
following minimax problem: 



In situations where one or more of the tuning parameters must be within user specified 
regions, the above minimax problem becomes a constrained minimax problem. 

Many commercially available software products can directly be used to find the poles of a 
transfer function and the roots of a polynomial, and to solve the minimax or constrained 
minimax problem as formulated above, such as (1) the function "pole" in the "Control 
System Toolbox", which can find the poles directly from the transfer function Q, (2) the 
function "roots" in MATALAB. which can directly find the roots of the characteristic 
equation b(z)=0, and (3) the functions "minimax" and "fminimax" in the "Optimization 
Toolbox", which directly provides solution to the constrained or unconstrained minimax 
problem formulated above, all of which are easy to use and commercially available from 
The Math Works Inc. 



min max{|zi|, |Z2|, .... |Zj|} 



